Non-Markovian master equations from entanglement with stationary unobserved 

degrees of freedom 
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We deduce a class of non-Markovian completely positive master equations which describe a system 
in a composite bipartite environment, consisting of a Markovian reservoir and additional stationary 
unobserved degrees of freedom that modulate the dissipative coupling. The entanglement-induced 
memory effects can persist for arbitrary long times and affect the relaxation to equilibrium, as well as 
induce corrections to the quantum-regression theorem. By considering the extra degrees of freedom 
as a discrete manifold of energy levels, strong non-exponential behavior can arise, as for example 
power law and stretched exponential decays. 
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I. INTRODUCTION 

Irreversible, dissipative quantum dynamics (such as of 
an open system embedded in an environment of uncon- 
trolled degrees of freedoms) differs drastically from re- 
versible dynamics described by a unitary time-evolution 

operator [l 0, II 0, IE IE 0]. An exaCt descri P tion of 
the dissipative dynamics can be given in the projector 
formalism 0], which results in a master equation for the 
reduced density matrix. In most cases, analytic progress 
can only be made under the Markovian hypothesis, which 
requires that correlations between the system and the 
environment decay faster than the characteristic inverse 
dissipation rate I/7. For weak coupling, Lindblad equa- 
tions can be derived which provide completely positive 
mapping of the density matrix from initial to final con- 
ditions; the underlying quantum dynamical semi-group 
can also be motivated from assumptions that are in- 
dependent of the weak-coupling requirement 0.0 A 
cornerstone of this framework is the quantum-regression 
theorem (QRT) 0,0, which relates multiple-time corre- 
lation functions to single-time expectation values. Fea- 
sible exact master equations beyond the Markovian ap- 
proximation are known when a spin S El or a har- 
monic oscillator |lfj is embedded in a bosonic bath. Also, 
transient deviations from Markovian behavior on times 
shorter than I/7 are well understood [Til IT2I Il3| . In gen- 
eral, however, only few general results and manageable 
models are known for non-Markovian dynamics beyond 
the transient regime [Tl ITl ITtl ITl 111 fis| . 

The few recent results about strong non-Markovian 
effects in quantum master equations were obtained in 
the context of complex environments [l5| . continuous 
time quantum random walks |16| . stochastic Hamiltoni- 
ans 01 > an d continuous measurement theory , as well 
as in more mathematical settings ^^|. In general, con- 
ditions for a QRT for non-Markovian dynamics have not 
been established in the past. 

In this paper, we present a new framework for the 
characterization of non-Markovian quantum system dy- 
namics. We show that long-time non-Markovian effects 



naturally occur in composite environments where the 
system-to-reservoir coupling strength depends on addi- 
tional quantum degrees of freedom. The non-Markovian 
behavior arises even when the reservoir itself can be de- 
scribed in a Markovian approximation because the me- 
diating degrees of freedom become entangled with the 
system degrees of freedom. These effects can persist for 
arbitrary long times, far beyond the transient regime. 
When the mediating degrees of freedom are eliminated, 
the dynamics can be rewritten as the statistical super- 
position of conventional master equations with random 
dissipation rate. For weak coupling we obtain a class 
of non-Markovian Lindblad equations which preserve the 
complete positivity of the solution map known from the 
Markovian case. However, we find that the QRT in gen- 
eral is not fulfilled, the only exception being the approach 
to a stationary limit which is independent of the random 
dissipation rate. We also present an effective approx- 
imation to the system dynamics, which facilitates the 
comparison with previous results for non-Markovian dy- 
namics. As an illustrative example we analyze the non- 
Markovian dispersive dynamics of a two level system. By 
assuming as the extra system a manifold of quantum lev- 
els, we demonstrate the possibility of getting strong non- 
exponential decays. 



II. COMPOSITE ENVIRONMENTS 

We start from a full microscopic description, consider- 
ing a system S that interact with a composite bipartite 
reservoir consisting in a bath B endowed with extra un- 
observed degrees of freedom U, which also participate in 
the system-environment interaction. The total Hamilto- 
nian reads 

H T = H S + (Hu + H B ) + XHj (1) 
with the tripartite interaction 

#/ = gs®(<2t/®Qs)- (2) 
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The identity XHj = (XQu) ® (qs ®Qb) implies that the 
operator Qu sets the system-reservoir coupling strength. 
For simplicity we assume that Qu is a constant of motion 



[Hu, Qi 



0. 



(3) 



Clearly, this assumption remains valid when the dynam- 
ics of Qu is slower than the dissipative relaxation 20] . 



Reformulation in terms of a random interaction 
parameter strength 

The dynamical evolution of the total density matrix 
p T (t) is given by 



p T (t) = exp[£ T £]pt(0), 



(4) 



where Ct[*] = (— i/K)[HT, •] is the total Liouville su- 
peroperator. In order to relate these dynamics to 
conventional dynamics in environments with fixed cou- 
pling strength, we eliminate the unobserved degrees 
of freedom Qu for usual factorizing initial conditions 
Pt(0) = ASB(0)<8>P;y(0) HU- The reduced density matrix 
PSB(t) = Tr;y[/9 T (i)] of the system S and the bath B is 
then given by 

Psb (t) = J2 Pr ex PK £ tf + c b + X R Ci)t] Psb(0), (5) 

R 

where Cr and Cr are the Liouville operators of the sys- 
tem and the bath respectively, and £/[•] = (— i/K)[{qs <8> 
Qb),*]- The index R runs over the eigenstates \R) of 
Hu- We introduced the probabilities 



P R = (R\ Pu (0)\R) 
and the weighted coupling strengths 
X R = X (R\ Qu \R) . 



(6) 



(7) 



Indeed, Eq. (J5J can be interpreted as a statistical average 
PSB(t) = (psb(^)) over solutions Pgg(t) corresponding 
to a Hamiltonian 



H' T = Hs 



H B + X R q s <8> Qb 



(8) 



with fixed interaction parameter A^. Each solution 
p§ B (t) participates with probability Pr. It follows that 
the reduced system density matrix ps(t) = Tr B [ps B (t)] 
can be obtained from the reduced density matrices 
Pn.(t) = ^bIpsb^)] by a similar average 



Ps{t)=J2 P nPR(t) = <P«(*))- 



(9) 



The random coupling formulation of the system 
dynamics allows to incorporate the previous knowledge 
about dissipative systems with fixed coupling strength 

QllHilliftiiHamHimp (it also may 
form the basis for efficient numerical simulations). For 



instance, if the map pr(0) — > pr(€) is completely positive 
then this is inherited by the map ps(0) — * Ps(t). 

In the rest of the paper, we will use the random de- 
scription to make further analytical progress for the case 
that the evolution of pR(t) can be described by a Marko- 
vian Lindblad equation p| . 



III. NON-MARKOVIAN LINDBLAD 
EQUATIONS 

When the correlation times of the bath B are the short- 
est time scale, a Markovian approximation applies. For 
factorizing initial conditions of the total density matrix, 
Asb(O) = ps(0) (8) Pb{0), and weak coupling of S and B 
(Xr <C 1), one then can describe the evolution of the re- 
duced density matrix pR(t) by a Lindblad equation 



dpRjt) 
dt 



C H [p R (t)} + lR C[pR(t)}. 



The random dissipation rate is given by 
1r = i(R\Q 2 u\R), 



(10) 



(11) 



where 7 is determined by the spectral density of the en- 
vironment evaluated at a characteristic frequency of the 
system. The Lindblad superoperator reads 

£M = JE([^'*^ + ^-•.^]), (12) 



with the operators {V a } acting on the Hubert space of S. 

From this description, it follows that the evolution 
of the reduced density matrix ps(t) = (p R (t)) is non- 
Markovian. In fact, note that due to the statistical cor- 
relation between 7^ and PR(t), the average of Eq. ijTU|) 
can not be written as a evolution which is local in time 



dpsjt) 
dt 



^ {C H + £'}[p s (t)}, 



(13) 



with £ being some extra Lindblad superoperator. Then, 
for obtaining the corresponding non-Markovian master 
equation, we first write the average Eq. @ in the Laplace 
domain as 



Ps(u) = 



1 



u - (C H + jCr) 



Ps(0) = {Gr(u)) P s(0), 



(14) 

where Cr = JrC and u is the Laplace variable. In order 
to cast this expression into a deterministic closed evolu- 
tion equation we have to interchange the average over the 
random dissipation rate and the operator-inverse in the 
definition of Gr. We employ the identity 



Ps(u) 



1 



(Gr(u){u-(C h +Cr)]) 



(G R (u)} Ps (0), (15) 
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and define a deterministic superoperator L by demanding 



Ps(u) 



Ps(0). 



(16) 



u - [C H + L(u)] 
The superoperator L is then determined by the condition 



{G R {u)C R ) = (G R (u))h(u) 
in the Laplace domain, and fulfills 



(17) 



(G R (t)C R )= dr {G R (t — t))L(t) (18) 
Jo 

in the time domain |22|. 

The consequence of this procedure is a deterministic, 
closed, non-Markovian evolution equation of the reduced 
density matrix, 

= C H [ Ps (t)]+J*drUt-T)\p s (T)]. (19) 

This equation has arbitrary long memory compared to 
the mean dissipation rate ("Jr). By construction from 
the average of Eq. (jTU|l . the solution map ps(0) — > ps(t) 
of Eq. iJTSj l is completely positive. 



A. Quantum-Regression Theorem 

For Markovian Lindblad equations the QRT 0, 0J pro- 
vides direct relations between expectation values of sys- 
tem observable and their correlation functions. We now 
use Eq. (|19fl to show that the theorem cannot be taken 
for granted for composite environments. Let us introduce 
a complete set of operators {A^} of the system, collected 
into a vector A, and consider the expectation values 



A(t) = TW[A(*)pr(0)], 
as well as the correlation functions 



(20) 



S(t)A(t + t)= Tr SU B[S(t)A(t + t) Pt (0)}, (21) 

where S is an arbitrary operator for the system. The 
time dependence of the operators refers to a Heisenberg 
representation with respect to the total Hamiltonian Q . 

Based on the random formulation (|9I1UI) of the dynam- 
ics, the expectation values and correlation functions can 
be written as averages over the dissipation rate, 



Affl - (Trs[Ap R (t)}) = (A(t) R ), (22a) 
S(t)A(t + T ) = (Tr s {Ae^+ c ^ r [p R (t)S]}) 

(22b) 



= (S(t)A(t + r) R ). 
The expressions deliver evolution equations 
d 



dt 



A(t) - (M R A(t) R ), 



(23a) 



r/r 



S(t)A(t + T) = {M R S{t)A(t + T) R ), (23b) 



where the matrix acts on the indices of A and is 
defined by the condition 

Tr s {A(£ H + Cr)[S}} = MrTt s {AS}. (24) 

When -fR is fixed, the evolution equations 123a|) for ex- 
pectation values and l|23bjl for correlation functions are 
identical, which recovers the QRT for Markovian dynam- 
ics. In the non-Markovian case, however, both equations 
still involve the average over the dissipation rate. 

By using the same procedure as for the density matrix, 
we can transform Eq. I|23a[) into a closed deterministic 
evolution equation, 



dt 



A(t) = - / drM(t - t')A(t'). 



(25) 



The deterministic kernel matrix M(i) fulfills an equation 
similar to Eq. I|17f) . but written in terms of and its 
corresponding propagator. 

Equation l|23b|) has the same structure as Eq. (|23a|l . 
but in the remaining average over the dissipation rate we 
are confronted with a subtlety. While Eq. (|23a|l is defined 
with initial conditions fixed at t = 0, Eq. (|23b|l gives the 
solution with initial condition S(t)A(t) R at finite time t. 
From the definitions in Eqs. (|22b|) and 124|) we find sta- 
tistical correlations between M# and S(t)A(t) R , which 
both depend on 7^. Dynamically these correlations can 
be understood by realizing that S(t)A(t) R is of the form 
of a single-time expectation value and hence evolves ac- 
cording to Eq. H23a[) ■ We still can disentangle the average 
over the dissipation rate by the procedure employed for 
the density matrix and the expectation value, but in- 
stead of a homogeneous equation of the form H25|) obtain 
an inhomogeneous equation 



dr 



S(t)A{t + t) 



dt' M(t - t')S{t)A{t + t') 



+I(*,r) 



(26) 



where the term I(i, t) accounts for the correlations. 

The QRT is fulfilled when the inhomogeneity I(t, t) 
vanishes, as is the case for Markovian dynamics where 
the average over the dissipation rate is absent. Equation 
(|22b(l implies that the inhomogeneity dies out in the long- 
time limit if the asymptotic state p R (oo) does not depend 
on 7a; the QRT is then asymptotically valid. However, 
if the asymptotic state pr(oo) depends on 7^ the inho- 
mogeneous term will contribute at all times, even in the 
asymptotic regime, and the QRT is invalidated. 



B. Effective Approximation 

In order to obtain a general characterization of the dy- 
namics, we introduce the following approximation. In 
Eq. (|17fl we discard the dependence introduced by the 
Lindblad superoperator C in the propagator Gr(u), i.e., 
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Cr — > — 7rI. Thus, we can write the approximated equa- 
tion as 



7i? 



which is solved by 



C 



U- Ch + 1R 



L(u) K(u-£ H )£, 



L(u), (27) 



(28) 



with the function 



7-R 



1 



(29) 



From here, the density matrix evolution reads 



dp s {t) 
dt 



C H \ps(t)]+ f dTK{t-T)e^ CH L[ Ps {T)]. 
Jo 

In this approximation all information about the extra 
system U is encoded in the kernel K(u), which is defined 
by Eq. (|29|l . While this approximation is not controlled, 
it is clearly useful for characterization of the possible non- 
Markovian effects. The solution of Eq. (|30|l will differ 
from the exact solution of Eq. (|I9(I only through small 
time dependent corrections, of order 1, of the decay rate 
parameters. On the other hand, we notice that the struc- 
ture of the evolution Eq. I|30[l is similar to that found in 
Rcf. [3 

in the context of a continuous measurement 

theory 



C. Stochastic state representation 

In the previous approximation, it is not clear whether 
the final evolution guarantees the completely positive 
condition. Here, by introducing a stochastic representa- 
tion of the dynamics, we proof that indeed this condition 
is preserved. 

In the Laplace domain Eq. I|30|l reads 

u Ps (u) - p s {0) = {C H + K(u - C H )C}p s {u). (31) 

By assuming that C — £—1 [23(, with £[•] = J2 a Va • V^J, 
we arrive at 



Ps{u) 



1 



1 - w(u - C H ) 



ps(0), 



(32) 



(33) 



I — w(u — Ch)£ J u — jCh 

where we have introduced the function 

i \ I Tfl 
w(u) = ( 

From here, it is possible to obtain the formal solution 

Ps (t) = P Q {t)e tc » Ps (0) (34) 
+ / dr w{t - T)e {t - T)CH 8[ Ps {t)]. 



where Pq(u) = [1 — w(u)]/u = ((u + 7_r)~ 1 ) . This equa- 
tion have a clear stochastic interpretation. It corresponds 
to an average of a stochastic density matrix p st (t) whose 
evolution consists in the application, at random times, 
of the superoperator £, while in the intermediate inter- 
vals the system state evolves with its unitary evolution 
exp[t£#]. The statistics of the random times is dictated 
by w(t), which can be interpreted as a waiting time dis- 
tribution, i.e., it is the probability density for an interval 
r between consecutive applications of £ . In correspon- 
dence, Po(t) = 1 — J Q dTiv(T) is the survival probability 
associated with w(t). Then, the first term in Eq. I|34|l 
represents realizations without any application of the su- 
peroperator £, while the integral term accounts for all 
other realizations. Thus, we can write ps{t) — ((p s t(t))), 
where ((•••)) denotes the average over the random times 
at which £ is applied. As each realization preserves the 
complete positivity, this property is also present in the 
averaged evolution. 

The previous stochastic framework allows us to clar- 
ify the role of the kernel K(t). This follows after intro- 
ducing the sprinkling distribution [24| f(t) = w(t)9(t) + 
JqW^ — t)/(t), where 9(t) is the step function [25| . 
From its definition, f(t) gives the probability density 
for an event at time t, disregarding the possibility of 
extra events in (0,£). In the Laplace domain it reads 
f(u) = w(u)/[l — w(u)]. From here it is simple to get the 
relation K(t) = df(t)/dt, which defines the kernel as the 
rate of the sprinkling distribution. 

When the unitary dynamics commutate with the ac- 
tion of the superoperator £, in an interaction represen- 
tation the stochastic dynamics reduces to that presented 
in Ref. ^(|. On the other hand, a similar stochastic 
interpretation may be proposed for the exact evolution 
Eq. I|19|) . involving many renewal processes. Neverthe- 
less, their specific form depends on the details of each 
problem. 

It is interesting to note that expressions similar to 
Eq. 134[) arise in the context of the micromaser theory 
|26l I27L |28| . This system consists in an electromagnetic 
cavity that is continuously pumped with excited atoms. 
In our scheme, the waiting time distribution w(t) can be 
associated with a non-Poissonian pump statistic, while 
the superoperator £ with the transformation produced 
in the cavity field by the passage of each atom. There- 
fore, the extra system U can be associated with the pump 
degrees of freedom. This comparison enlightens the dy- 
namical origin of the non-Mar kovian effects. 



D. Dephasing of a two- level system 

As an illustrative example of our results we con- 
sider a two-level system, described by the Hamiltonian 
Hs = (l/2)hujAO'z, where a z is the Pauli z-matrix, 
which is weakly connected to a composite environment. 
The reservoir B is described by the dispersive Lindblad 
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operator 



»cr 2 ]), 



(35) 



and the mediating part U is described by an arbitrary 
set {7_r, Pr} of dissipation rates and weights. 

The evolution of the density matrix is given by the 
non-Markovian Lindblad equation l|19|) . From Eq. 1171) 
we find the superoperator L(u) = K{u — Ch)£, with the 
kernel K (u) defined in Eq. (|29() . Thus, in this case, the 
evolution Eq. l|5U|) is exact. 

The density matrix solution can be easily found in 
an interaction representation with respect to the system 
Hamiltonian. We find the completely positive map 

Ps (t) = g+(Jt) P8 (p) + g-(t)a zPs (0)(T z , (36a) 



g ± (t) = -[l±Po(*)]. 



(36b) 



The function Po(t) = Pr exp[— ^Rt] is the survival 
probability defined previously. Depending on the dis- 
tribution of the dissipation rate, arbitrary forms of the 
decay can be obtained from this average over exponen- 
tial functions. Hence the non-Markovian behavior can be 
observed in the relaxation of the density matrix to the 
stationary state. 

Let us now illustrate the consequences for the QRT, 
hence, the expectation values and correlators of the vec- 
tor of Pauli operators A = {a x , a y ,a Z} I}. For the expec- 
tation value (in the interaction representation) we find 



A(t) = G(t)A(0), 

where the matrix propagator reads 
diag{Po(t), Po(t), 1, 1}. Consistently with 
the correlation functions 



G(i) 
Eq. 



S{t)A(t + t) = G(r)S(t) A(i) + l h(t, t). 



(37) 



(38) 



feature an extra inhomogeneity, which is given by Io = 
BTr s {[/3 S (0) - p 5 (oo)]SA}, with the matrix D = 
diag{l, 1,0,0}, and the function h(t,r) — Po(t + r) — 
Po(t)Po(T~). In the Markovian case h vanishes since the 
survival probability Pq (t) is an exponential, and the QRT 
is valid for all times. For non-Markovian dynamics, the 
theorem is valid in the long-time asymptotic, since the 
inhomogeneous term dies out as the equilibrium state 
ps (00) is attained [consistently, in this example pr(oo) 
is independent of the dissipation rate 7_r] . 



IV. ENTANGLEMENT WITH A DISCRETE 
MANIFOLD OF ENERGY LEVELS 

Up to now we have left unspecified the unobserved de- 
gree of freedom U, which determine the set {7_r, Pr}. In 
this section we will analyze the case in which U is defined 
by a discrete set of energy levels. 



In the effective approximation, the properties of the 
unobserved degree of freedom are introduced through the 
kernel K(t), which is associated with the renewal process 
defined by w(t) = ^2 R Pr ~iRe~~ lRt . This process has two 
characteristic times scales 



(39) 



When these constant are well defined, the average rate 
(7) defines an exponential decay of the waiting time 
distribution at short times [('y)t < 1], lim tl ^ 00 it;(w) ss 
{y}/(u + (7)), and the average waiting time (r) = 
J dttw(t) characterizes an exponential decay of w(t) in 
a long time regime [t > (r)], lim u ^ow(u) « 1 — u(t). On 
the other hand, in terms of the sprinkling distribution we 
get ]im t _ > o+/(t) = (7} and limt-oc/C*) = 1 /( T )- 



A. Entanglement with a two-state system 

First we assume that the unobserved system U is a 
two level system. The waiting time distribution then 
reads w(t) = Pt 7|e _7Tt +Pi 7|e _7i *, with the condition 
Pf + P[ = 1 and arbitrary rates 7f/j.. Introducing the 
rates rj = P T7i + P i7T and f3 = [{j 2 ) - {l) 2 ]/{l), the 
Laplace transform w(u) can be written as 



w(u) 



(7) 



u + (7) + (3<t(u) ' 



(40) 



with cr(u) = u/[u+i-]/(('j)(t})]. The corresponding kernel 
Eq. results as 



K(u) 



(7) 



(41) 



which in the time domain reads K(t) = {j)[S(t) — f3e~ nt ]. 
We note that the fluctuation rate (3 controls departure 
from a Markov kernel. The sprinkling distribution results 
as f(t) = (j)6(t) - [(7) - (t)- 1 ]^ - e'^), where we 
have used r)[l — ((7)(t)) -1 ] = (3. As expected, (7) and 
(t) _1 give respectively the asymptotic values of f(t) in 
the short and long time regimes, while rj gives the rate 
for the transition between these two regimes. 



B. Entanglement with a N-manifold of states 

Now we characterize the case in which the system U 
consists in a manifold of N states [0 < R < N — 1] 
whose consecutive energy difference is constant. We as- 
sume that the coupling strength of each level with the 
system-bath-set decreases in an exponential way with the 
level energies, as well as their stationary populations 



-1r = 7exp[-6P], Pr 



7T-^exp[-aP]. (42) 
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The constant 7 characterizes the Markovian decay of 
the system, and b is a free dimensionless parame- 
ter that measure the strength of the coupling be- 
tween each states and the system-bath set, i.e., Qu = 
Y^rZo ex P[~ bR/2]\R)(R\. The second free parameter a 
measure the exponential decay of the populations. By 
taking a — fttuo/kT, where fkoo is the difference of energy 
between consecutive levels, the populations correspond 
to a thermal distribution at temperature T. 
The average rate reads 



{7} = 7 



1 -e~ a 
1 - e -( a + b ) 



1 _ e -(a+b)N 



1 - e- a 



N 



(43) 



from where the average waiting time follows immediately 

1 — e~ a \ (\ — e -( a - b ) N 1 



[ ' 1 \l~e-( a - b ) J \ l-e- aN ) 



(44) 



In an intermediate regime, (7) -1 < t < (t} _1 , the 
corresponding waiting time distribution may present a 
power law behavior. In fact, in the limit of N — > 00, 
where the TV-manifold states is equivalent to a thermal 
harmonic oscillator in equilibrium at temperature T, it 
is possible to proof that, after a time transient of order 
I/7, the waiting time distribution behaves as |29| | 



w(t) « l/( 7 t) 



l+Q 



a/b. 



(45) 



For finite N, this behavior is also present, nevertheless 
the asymptotic behavior changes to an exponential decay 
with rate (r} _1 . 

For a < b [0 < a < 1] the behavior of w(t) can be 
captured with a simple analytical expression. Taking in 
account the results for the case N = 2, Eqs. (liT) )) - (|4*T1l . 
we propose the complete monotone function |16( 



w(u) 



(7) 



+ ( 7 }+/3 1 -«a a (w) : 



(46) 



where as before (3 w [(7 2 ) — (7) 2 ]/(7), and now a a (u) — 
[(w + 7 c ) Q —7"]- Using the relation K{u) — w(u)/P n (u), 
the kernel reads 



K(u) 



<7) 



1 + /3^ a a a (u)/u 



(47) 



From lim tl _ + o'u ; ( u ) ss 1 — [u + f3 1 ~ a o- a (u)]/('~f), the asymp- 
totic exponential decay of w(t) with rate (t) -1 can be 
fitted with the cutoff j c after imposing the relation 
a (/3/7c) l Q = (l)( T )~ 1- On the other hand the constant 
P fits the power law regime. In fact, in the limit N — > 00, 
the average waiting time Eq. (|44|l is infinite, which im- 
plies j c = 0. Thus, lim u ^ w{u) « 1 - /3 1 - a w^/(7), 
implying a pure power law asymptotic behavior |3Cj ]. 

The previous analysis demonstrates that the transient 
behavior between the short and asymptotic exponential 
decays of w(t) is described by a power law. As we have 
seen in the previous section, this behavior is in general 
reflected by the system dynamics. On the other hand, we 



notice that when N — > 00, by maintaining the average 
rate (7} fixed, the fluctuation rate (3 reaches it maximum 
value in the limit of both, strong coupling 6 « 1 and 
small populations decay a«l (which can be also read as 
a high temperature limit). In this case, after the transient 
t < (7) , the waiting time distribution and the kernel 
can be approximated by the expressions 



w(u) 



A, 



K(u) » A c 



(48) 



with A a w (7} //3 1_Q w (7)"- These expressions cor- 
respond to a fractional derivative evolution |l6l |. where 
stretched exponential and power law behavior arises 
jointly. 



V. CONCLUSIONS 

We have shown that a system in a composite bi- 
partite environment (where the system-to-reservoir cou- 
pling depends on other degrees of freedom) follows non- 
Markovian dynamics even when the reservoir itself can 
be eliminated by a Markovian approximation. The non- 
Markovian effects originate in the entanglement of the 
system with the mediating degrees of freedom, and may 
persist for arbitrary long times. 

Our results are derived from a random rate reformula- 
tion of the dynamics in the composite environment which 
allows to make full contact to the established theory 
of dissipative systems with constant coupling. On this 
basis, we formulated non-Markovian Lindblad equations 
which provide complete positive mappings of the density 
matrix from initial to final conditions, and identified con- 
ditions for the quantum-regression theorem. It should be 
noted that the random rate formulation is not restricted 
to the Lindblad master equations but can be applied to 
other master equations, including exact master equations 
which may already include non-Markovian effects at fixed 
coupling. 

In an effective approximation all information about 
the extra degrees of freedom is introduced by a mem- 
ory kernel. The corresponding density matrix evolution 
can be interpreted in terms of stochastic process in the 
system Hilbcrt space. System decay behaviors ranging 
from stretched exponential to power law can be obtained 
by taking the system U as a discrete manifold of states. 

The present formalism substantiat es p revious results 
on non-Markovian master equations |14L IT5L ITtj ITtI ITsl 
Il9| , and puts them into the alternative and greater per- 
spective of systems embedded in a composite environ- 
ment. Our motivation to study this kind of environment 
arises from recent ex per iments on fluorescent single quan- 
tum dots mHHHll, where non-Markovian effects on 
time scales much larger than 1/7 were found. While the 
underlying physical mechanisms are not completely clear, 
it has been argued [35l |3(j that the experimental results 
can only be recovered when one accounts for additional 
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degrees of freedom which modulate the dissipative cou- 
pling. Thus, besides its theoretical interest in the context 
of strong non-Markovian effects in open quantum sys- 
tems, the characterization of the dynamics induced by 
composite environments may be also of interest in those 
experiments. On the other hand, we believe that the 
present results may be useful in modeling the dynamics of 
open quantum systems embedded in complex structured 
host environments. In fact, the tri-partite interaction in- 
vestigated by us naturally arises when one considers the 
dynamical effects of a disordered condensed-matter en- 
vironment on a system coupled to a (say, Markovian) 
environment. The coupling strength of the system to the 
Markovian environment is proportional to the density of 
states of the Markovian environment at a characteristic 



frequency of the system; this density of states, in turn, 
depends on the condensed-matter environment (e.g., the 
disorder configuration, or the charging of trap states). 
Taking the dynamics of the condensed-matter environ- 
ment into account, the coupling strength then becomes 
as a dynamical variable, leading to the multiplicative tri- 
partite coupling that was considered in this paper. 

In the description of condensed-matter systems, mul- 
tipartite interactions also frequently appear as a conse- 
quence of a suitable transformation (such as the Frohlich 
transformation) |37| . Such transformations can be used 
to eliminate (at least to some order of a small parame- 
ter) the time-dependence of a certain subsystem, which 
then takes on the characteristics of our (stationary) un- 
observed degrees of freedom. 
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